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Abstract When a nonequilibrium growing interface in the presence of a wall is considerad a nonequi- 
librium wetting transition may takc place. This transition can be studied trough Langevin equations 
or discrete growtli models. In the first case, the Kardar-Parisi-Zhang equation, which defines a very 
robust universality class for nonequilibrium moving interfaces, with a soft-wall potential is considered. 
While in the second, microscopic models, in the corresponding universality class. with evaporation and 
deposition of particles in the presence of hard-wall are studied. Equilibrium wetting is related to a 
particular case of the problem, it corresponds to the Edwards-Wilkinson equation with a potential in 
the continuum approach or to the fulfiUment of detallad balance in the microscopic models. In this 
review we present the analytical and numerical methods used to investígate the problem and the very 
rich behavior that is observed with them. 
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1 Introduction 



Wctting [1] is well cxemplified by considering a liquid droplet on a substrate. Depending on the physical 
properties of the system, determining the shape of the droplet, the substrate will be more or less wet. 
More specifically, the contact angle & (see Fig. [T|) is related to the surface tensions trough Young's 
equation, 

eos O = {(Tsy - crs,L)/c^L,v, (1) 

whcrc (Ts.Vj o's.l and ctl.v sltc the surface tensions of the substrate- vapor, substratc-hquid and Uquid- 
vapor surfaces, respectively. Total wetting happcns if 6* = O, and O < < tt corresponds to partial 
wetting. 

Cahn [2] observed that by approaching the critical temperature Te, for T < Te, the liquid- vapor 
surface tensión aL,v goes to zero faster than the difference (erg y — gs^l), therefore, at a temperature 
Tw < Te a wetting transition should take place. In this wetting transition, 6* > O for T < Tw and 
= for T > Tw- This observation introduced the notion that wetting could be viewed as a type of 
critical phenomena. 

More generally, a wetting transition occurs in a thcrmodynamic system constituted of a bulk phase 
A and a substrate that attracts a second coexisting phase B. Depending on the control parameters, 
like e.g. temperature and chemical potcntial, the system will be in the moving phase or in the bound 
phase. At the wetting transition a macroscopic layer of the absorbed phase B is formed and phases A 
and B coexist. In the bound phase only the phase A is stable and the AB interface stay pinncd to the 
wall, whereas in the moving phase the AB interface grows and only phase B is stable. 

Abraham [3] introduced a two-dimcnsional Ising modcl that could be solved exactly and displayed 
the same kind of transition prcdictcd by Cahn [2]; it was also shown that in the solid on solid (SOS) limit 
the model still cxhibitcd a wetting transition. The advantage of taking this limit is that it simplifies the 
calculations and with it scvcral exact results can be obtained using the transfcr matrix method [4-6] . 

The SOS model, defined on a one-dimensional lattice with pcriodic boundary conditions and size 
L, has the foUowing Hamiltonian, 

L L 

H = jY,\K+^-h,\+Y,V{h{), (2) 

i=l i=l 

where /i¿ > O is a discrete random variable representing height, J is a positive coupling, V{hi) is 
a potential accounting for the interaction between the wall = 0) and the absorbed phase and 
no overhangs are allowcd. The región x > hi corresponds to the phase A and the región x < /i¿ to 
the coexisting phase B, therefore, the height gives the position of the AB interface. By varying the 
temperature this modcl presents a wetting transition at a temperature Tw, at which the average height 
of the interface diverges. 

On a coarse-grained level equilibrium wetting transitions can be studicd using the Hamiltonian [7] 

H = j dMf|VMx)P + ^[Mx)]}, (3) 

where x is a continuous variable giving the position in the d-dimcnsional substrate. The height /i(x) 
gives the position of the interface and the potcntial V [ft-(x)] accounts for the presence of the substrate 
and a possible interaction between it and the interface. In ordcr to consider the dynamics of wetting, 
Lipowsky [8] introduced the foUowing Langcvin equation, 

^ - - + -V^Mx, t) - + C(x, t), (4) 

where the deterministic part of it originates from —5H[h{x)]/&h{x), with H given by and C(x, í) 
is a Gaussian noise with zero mean and variance given by 

(C(x, í)C(x', í')) = D5{^ ~ ^r-'S{t - t'). (5) 
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Fig. 1 Two drops on a substrate. The left one has a contact angle > 7r/2 and wets less the substrate, in 
comparison to the right one with O < n/2. 



Equation (¡4|) is the Edwards-Wilkinson equation [9], which is known to describes the motioii of an 
equüibrium interface with vclocity a (see [10,11]), with an extra term: the potential í/[ft,(x)]. 

The Kardar-Parisi-Zahng (KPZ) [12] equation differs from the EW equation by the presence of 
the nonhnear term X{\7h)^. The basic difference of an interface described by the KPZ equation (when 
compared to one described by the EW equation) is that the velocity of the interface depends on the 
local interface slope V/i(x). The non-hnear term comes from an expansión of the velocity with respect 
to the local interface slope, and it can be shown that higher order terms in the expansión are irrelevant 
under a renormalization group transformation [10,11]. That is why the KPZ equation defines a robust 
universality class for noncquilibriuni nioving interfaccs. 

Since the paper of Cahn [2], equihbrium wetting has been intensively studied, theoreticaUy and 
experimentally (see [7,13,14] for reviews). Differently. nonequilibrium wetting is in it's very beginning 
and experimental studies on it still lacks. One way to study it is to add the nonhnear term of the KPZ 
equation A(V/i)^ to equation (¡4]), such equation was considered for the first time by Tu et aL [15]. 
Hinrichsen et aL [16] and Muñoz and Hwa [17] introduced microscopic models in the KPZ universahty 
class in the presence of a hard wall, showing that nonequilibrium wetting could also be studied with 
them. A considerable amount of work on nonequilibrium wetting has already emerged and our aim is 
to review it here. 

In this article we present a detailed description of nonequilibrium wetting, reviewing the very rich 
phenomenology that is obtained from numerical and analytical methods used to study the Langevin 
equation and the microscopic models. Two reviews on the subject, more focused on the Langevin equa- 
tion approach, are [18,19], and there is some overlap between them and the present one. Nevertheless, 
in the account wc make here many aspects of nonequilibrium wetting are prescnted in more detail, 
mainly when it comes to microscopic models. We also obtain a few new results that will be pointed 
out when they come. 

The review is organized in the foUowing way. In the next section we define the Langevin equations 
and the microscopic models to be analyzed. Moreovcr, the quantities of interest and the exponents 
associated with them are defined. Sec. 3 contains exact calculations, that are possible for microscopic 
models when detailed balance is fulfiUed. In Sec. 4 Monte Cario simulations of microscopic models 
and numerical integration of the Langevin equations are presented. The content of Sec. 5 is mean 
field approximations for microscopic models. The case of a substrate with more than one dimensión 
is treated in Sec. 6, where mean ficld approximations for the continuum model, naive power-counting 
and some renormalization group argumcnts are used. We end with final remarks in Sec. 7. 



2 Definition of the problem 



2.1 Contimmm model 



The Langevin equation we are going to consider, describing the niotion of an interface in the presence 
of a wall, is the bounded KPZ (bKPZ) equation, givcn by [15] 

= a - ¿Vih) + ^V2/i(x, t) + A(V/i(x, í))2 + C(x, t). (6) 
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This is cquation ^ added with thc nonlincar tcrm A(V/i(x, í))^. In contrast to thc unbound KPZ 
interface, the sign of A is of great importance, because, as we will see, it leads to different univcrsality 
classes. We cali them bKPZ+ and bKPZ- universality classes. 

The potential V{h) has to account for the presence of a substrate and the interaction between the 
absorbed phase and the substrate. A natural way to consider the presence of a substrate is to forbid 
negative heights with a hard-wall potential. that is givcn by 



The problem with this potential is that it is not suitable for calculations (analytical and numcrical). 
A potential that overcomes this problem and is used in thc study of cquilibrium and nonequilibrium 
wetting with short range interactions is [7] 

V{h) = - cxp(-sh) + ^ cxp{-2sh), (8) 
s 2s 

wherc s > 1 controls the hardness of the wall, 6 is a control parameter that in equilibrium and near the 
transition is proportional to \T — Tw\ and c > 0. If 6 > O the wall is repulsive and the term exp(— 2s/i) 
is irrelevant, with 6 < O and c > O we have an attractive wall. The order parameter of the wetting 
transition is n = e~^. In the bound phase thc mean hcight is finito and thc interface stays pinncd to 
thc substrate, therefore n > 0. In thc moving phase thc mean hcight grows linearly with time and 
n — 0. As we will see, in the microscopio models a hard-wall is considcrcd, and thc critical behavior of 
thcm is thc same as thc critical behavior of thc bKPZ cquation with thc soft-wall potential ([8]). 

It turns out that thc bKPZ cquation (¡S]), with V(h) givcn by (¡S]), can be transformcd into a Langevin 
cquation with multiplicativo noisc. This is done with the Cole-Hopf transformation n = e^'\ which, 
in the case that thc sign of thc nonlincar tcrm is negative, gives the following cquation for the order 
parameter [20], 

^^^^ = -¿Vi^) + í^V^níx, t) + n(x, í)C(x, t), (9) 

wherc without loss of gcncrality we set A — —a and chosc to intcrprct thc Langevin cquation in the 
Stratonovich scnsc (intcrpreting it in the Ito sense would just produce a shift in the factor multiplying 
thc linear term [21]). The potential, as a function of n, is now givcn by 

V{n) = -n^ + -^n^+' + ^—n^+^', (10) 
^ ' 2 2 + s 2 + 2s ' ^ ^ 

wherc the linear term —an appearing in cquation ([9]) is now incorporated in the potential. This cquation 
was introduccd by Grinstein et al. [20], latcr it was pointcd out that it was cquivalcnt to thc bKPZ 
cquation [15]. The universality class it defines is callcd multiplicative noise 1 (MNl), obviously MNl 
and bKPZ- are thc same universality classes. Thc MNl cquation is also rclated to the synchronization 
transition in couplcd map latticcs [22], therefore, there is a relation between this transition and the 
wetting transition. Wc do not discuss this point here, references about it are [18,23,24]. 

By applying thc same transformation to thc bKPZ+ case (now with A = cr), thc rcsulting cquation 

is 

dn(x,t) d ^. . „9 , , „ (Vn(x,í))-^ / n >, ^ / x 

l,; ' = -:r^in) + aV^n X, t) - 2a^ ^^-^ + ?7, x, í C x, í . 11 

Ot dn n(x, í) 

This cquation is known as thc multiplicative noisc 2 (MN2) cquation. The MNl ^ and MN2 (fTT|) 
equations will be very important in thc course of this review, because they are more convenient in 
several situations. For cxamplc, integrating them numerically is simplcr than intcgrating thc bKPZ 
cquation dircctly, which suffers by numcrical instability problcms [25]. 

With this map between the bKPZ equations and the MN equations it is casy to see what are 
thc cffccts of considering an uppcr wall instead of a lower wall. An upper wall is implemented with 
thc potential V{h) = ^exp{sh) + j^exp{2sh) in cquation ([6]), in this case an negative (positivo) A 
corresponds to the case of an lower wall with positivo (negative) A. This can be verified by pcrforming 
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Fig. 2 Transition rates for RSOSW model (left) and example of an interface configuration with transitions 
that may take place. 



the Cole-Hopf transformation, but now with n = e'', in the bKPZ equation with an upper waU. The 
result is: with a positive (negativo) A the obtained equation is MNl (MN2). In this review we always 
consider a lower wall, also in the microscopic models. 

Bcforc going to the definition of the microscopic models we point out that, although the bKPZ 
equation is defined for d spatial dimcnsions, in most of this review we will restrict to the case d = 1, 
the exception being Sec. 6. For this reason, the microscopic models are defined, in the foUowing, on a 
one-dimensional substrate. 



2.2 Microscopic models 



2.2.1 Restricted solid on solid model 



The model for nonequilibrium wetting described in what foUows was introduced by Hinrichsen et 
al. [16]. It is a growth proccss taking place in a one-dimensional discrete latticc with periodic boundary 
conditions and size L. To each site i a random variable hi is attached, it can take the valúes hi = 
0, 1,2,3... and is interpreted as the interface height. Nearest ncighbors respect the restricted solid on 
solid (RSOS) constraint, i. e., 

\h,-h,±i\<l, (12) 

which introduces an effective surface tensión. The interface evolves in time by random-sequential up- 
dates in the foUowing way. A site i oí the lattice is randomly chosen and the processes that may occur 
are (see Fig. [2]): 

(a) deposition of a particle (/i¿ — > hi + 1) with rate q, 

(b) evaporation of a particle {hi hi — 1) at the cdges of platcaus with rate r, 

(c) evaporation of a particle {hi ^ hi — í) from the middle of a platean with rate p. 



If the final configuration would viólate the RSOS condition or would lead to a negative height, than 
it is not carried out. After L attempts a Monte Cario step is completed, and time is increased by one. 
The initial condition is a flat interface at height zero and we, without loosing generality, set r = 1. A 
hard wall is prcsent because evaporation events at the bottom layer {hi = 0) are forbidden. Since this 
model respects the RSOS condition and thcrc is a hard-wall, we ñame it RSOSW model. 

The wetting transition can be explained as foUows. Consider a free interface (negative heights are 
allowed). It may propágate, depending on the deposition and evaporation rates, in the directions of 
increasing or decreasing height. In the phase where the velocity of the interface v is positive (increasing 
height direction), the presence of the wall makes no difference: after some transient the interface will 
propágate with the same velocity as if the wall was not present. In contrast, when the rates are such 
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Fig. 3 Phase diagram of the RSOSW model. At p = O the transition is in the direct percolation (DP) 
universality class; for O < p < 1 it is in the bKPZ- universality class; for p = 1, where detailed balance is 
satisfied, it is in the bEW universality class; for p > 1 in the bKPZ+ universality class. 



that w < O, the wall changes the scenario completely because the interface stays bounded to the waU. 
Therefore, by forbidding negative heights, a wetting transition from a bound to a moving phase takes 
place, and the phase transition hnc qdp) corresponds to a frcc interface with w = 0. The phase diagram 
of the RSOSW model is displayed in Fig. H 

This model was generalizcd by Hinrichsen et al. [26], in order to include an attractive interaction 
between the substrate and the interface. This is done by considering a different deposition rate < q &t 
zero height. Because is smaller than q the detachment of the interface from the substrate becomes 
harder and, therefore, this change in the dynamical rules simulates an attractive forcé between the 
substrate and the interface. As we will show, using microscopic models and the bKPZ equation, the 
presence of an attractive potential Icads to new physics. 

The sign of the factor multiplying the nonlinear term of the the KPZ equation can be determined in 
a microscopic model. Considering a initially tilted interface, it is related to how the interface velocity 
varies with the tilt [10]. In Fig. [3] the doted-line where A = O is displayed, above (below) it A < O 
(A > 0). This shows that for O < p < 1 the phase transition is in the bKPZ-, for p = 1 in the bEW and 
for p > 1 in the bKPZ+ miiversality class. p = O is a particular case that will be addressed in Sec. 4. 

2.2.2 Single-step model 

Another microscopic realization of nonequilibrium wetting is the so-called single-step model with a wall 
(SSW) introduced by Ginelli et al. [27]. Here we use a more specific versión of it, studied in [28,29]. 

In the SSW model, the difference of height of two neighbors sites is restrictcd through the condition 



As in the previous model, it is defined on a one-dimensional lattice of size L and with periodic boundary 
conditions. It evolves random-sequentially by the foUowing rules: 

(a) deposition of a particle (/i¿ hi + 2) with probability p, 

(b) evaporation of a particle {hi s- hi — 2) with probability 1 — p, 

at a chosen site i occur if the configuration after the move is carried out does not viólate the single-step 
constraint. The initial condition is hi — 1 {hi = 0) ií i is odd (even). In the case of the SSW model A 



\K - h,+i\ = 1. 



(13) 
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Fig. 4 Example of an interface configuration of the SSW model with some transitions that can occur. The 
horizontal hne indicates the height of the moving wall h. On the right is the interface after the substrate moves, 
h^h+1. 



is proportional to 1/2 — p, thcrcforc for p = 1/2 we have the bEW case and for p > 1/2 {p < 1/2) the 
bKPZ- (bKPZ+) case. The incorporation of a wall, in comparison to the previous case, is a bit more 
complicatcd, it is done in the foUowing way. Different from the RSOSW model, it is not possible to 
vary the velocity of the interface keeping A fixed in the SSW model because there is just one control 
parameter, namely p. Fortunately, the velocity of the interface, in the long time limit, for the single-step 
model (free interface case) is known exactly; it is given by 



Therefore, in order to study nonequilibrium wetting the systcm is tuned to criticality by considering a 
wall that moves with velocity vl, given by (jl4|) . This is done by forbidding evaporation events below 
the substrate height h{t). 

In this paper we attain ourselves to the cases p ~ 1/2, p = 1 and p = Q which are in the bEW, 
bKPZ- and bKPZ+ universality classes, respectivcly. For p = 1/2, sincc w¿ = O, the height of the wall is 
h{t) = 0. At p = 1, during a simulation, after every At — 2(1 — 1/L) we increment the substrate height 
by one unity, /i — > h + 1. This means that all sites below the new h have their height increased and 
evaporation events at the new h are forbidden. For the bKPZ+ case {p ~ 0) after every Z\í ~ 2(1 — 1 /L) 
we have h h — 1. 

The SSW model is particular useful to obtain critical exponents numerically, since in this case 
one knows the critical point exactly. On the other hand the advantage of the RSOSW model is that it 
presents a richer behavior. The SSW model is also important in the study of mean field approximations 
(see Sec. [5]). 



2.3 Observables and related exponents 

Finally. we define the physical quantitics, that we are going to consider in this review, and the exponents 
associated with them. They are defined in the context of the RSOSW model and will be generalized 
to the continuum model at the end of this section. The generalization to any other microscopic model 
is straightforward. 

2.3.1 Scaling exponents 

A key quantity in what foUows is the mean height of the interface, defined by 



where () means the ensemble average. Another observable of great importance is the interface width. 



VL = (p-l/2)(l + l/L). 



(14) 



L 




(15) 



given by 



L 




(16) 
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Fig. 5 Interface width w, for the RSOSW model, as a function of time in the bound phase (left) and in the 
moving phase (right) for L = 64, 128, 256. In the bound phase q = 0.9 and in the moving phase q = 3.0, where 
p = 1.0. One can see that the saturation valué of w tends to a constant in the bound phase (smooth interface) 
and grows with L in the moving phase (rough interface). 



Considering an infinite system, if the interface width grows with time and does not reach a stationary 
valué in the long time limit, the interface is called rough. Otherwise, if w saturates after some transient 
it is smooth. In the one-dimensional case, the wetting transition is a roughening transition, the interface 
is smooth for q < Qc and rough for q > q^. Obviously, if the system is ñnite, even when the interface is 
rough, w saturates. When one wants to verify if an interface is rough or smooth, by doing simulations 
in finite systems, one has to consider different sizes and see how the saturation valué Ws{L) varies with 
L. If Wa{L) grows with L, than the interface is rough; if it tends to a saturation valué independent of 
L, then the interface is smooth. This is shown in Fig. [5l 

The KPZ universality class describes the self-afñne properties of the roughening interface under 
scale transformations. The scaling exponents are defined by the relations [10,11] 

w{t)^t\ í. (17) 

where tg is the time at which the interface width saturates. 7 is known as the growth exponent, a 
is the roughness exponent and z the dynamical exponent. They are not all independent, but rclated 
because of the Family- Vicsek scaling relation [30] 

w{t,L) = L'^f{tL-'), (18) 

where ,f{x) is a scaling function. From this last relation foUows that z ~ a/f3. In one dimensión the 
scaling exponents of the EW universality class are a = 1/2, /3 = 1/4 and z = 2; while for the KPZ 
universality class they are a ~ 1/2, (3—1/3 and z — 3/2 [10, 11]. 

We stress that the definition of the scaling exponents does not depend on the presence of the wall, 
they are rclated to the invariance of scale of w. However, the definition of critical exponents just make 
sense with a wall. As we will show, by introducing a wall, only one new independent critical exponent 
arises, the others can be determined by scaling relations and the valúes of the scaling exponents. 

2.3.2 Critical exponents 

The order parameter of the wetting transition is the density of sites at zero height, i. e., 

L 

Po = (¿"'^áh.,o). (19) 

i=l 

In the bound phase po > O, while in the moving phase, where the interface detachcs from the wall, 
Po = 0. Another order parameter for the transition is the velocity of the interface, which is zero at the 
bound phase and non-zero at the moving phase. 
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Fig. 6 On the left po{t) below {q — 0.97), at {q = 1.0) and above {q = 1.01) criticality, with L — 512 and 
p = 1.0. Note that at the critical point po{t) decays algebraically, above it goes to zero faster exponentially 
and below it saturates. On the right, typical bounded and moving interfaces are showed for difFerent times. For 
q = 0.97 is stays bounded to the wall and for q = 1.01 it advances. 



Fig. [ü] shows the typical time evolution of po above, below and at criticality. For q < Qc the order 
parameter reaches a constant positive valué in the stationary state. At the critical point, it goes to 
zero with a power-law behavior and for q > Qc it vanishes exponentially. The exponent 9 is defincd at 
the critical linc by the rclation 

po(í) ^ t-'. (20) 

Near and below criticality we have 

Pt>-{qc-qf, (21) 

whcrc Pq is the saturation valué of the order parameter. As discussed above the mean hcight is finito 
in the bound phasc and diverges at criticality. The exponent associated with this divcrgcncc is defincd 

by 

{h)^{q,,-q)-<, (22) 

where the above relation is valid near and below criticality. In one dimensión the interface width is also 
finite for q < q^ and diverges at the critical point. Since it has the same dimensión as {h) we expect it 
to diverge with the same exponent C^. In this scnse, whcn dealing with a onc-dimcnsional substrato, we 
consider the above definition also valid for w. When analytical calculations are possible, we calcúlate 
both quantities and show that, indeed, they have the same critical behavior. 

At the wetting transition, like the interface width, the mean height foUows the relation (/i) ^ . 
In the moving phasc {h) grows lincarly with time. Therefore, the interface vclocity goes to zero as the 
critical point is approachcd from above. The critical exponent /?„ is defincd by 

V^{q-qc)^\ (23) 

where the above relation is valid above and near criticality. 

The spatial correlation length diverges near the critical point as, 

í^--{qc-q)-"^. (24) 
The same happens with the temporal correlation length ^y, 

Cii ~ (gc-g)-"". (25) 

The saturation of the interface width w in a finite system, for a rough interface, happens when 
becomes of the same order of the system size L, at the time tg proportional to . AU this lead to 

íll-ei, (26) 
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which gives the foUowing scaling relation, 



Z = í^||/:/_L. 



(27) 



Another scaling relation is 



(28) 



it comes from (PÜ¡) . ([?!]) and (|25p . Since the spatial correlation length is of the order of the system size 
when w saturates, Wg ~ i" can be written as Wg ~ C"- Henee, with ((22)l and ([24]), we have 



We are considering the case qo = q where there is no attraction bctween the wall and the absorbed 
particles. For the bKPZ equation, where the order parameter is n = e^'', we have a non-attractive wall 
if b assumes some positive fixed valué (in this case the term exp(— 2s/i) is irrelevant and one can sct 
c = 0). The control parameter a is analogous to q in the RSOSW model and, therefore, relations (j20p 
and (|2ip becomc: n{t) ^ (valid at a = üc) and rig ^ (a^ — a)^ , where is the saturation valué. 

2.3.3 Attractive substrate 

If there is an attractive forcé between the substrate and the particles, different physical properties 
are observed. Within the RSOSW model the situation is as foUows. For the bEW and the bKPZ+ 
universality classes (p > 1), as the deposition rate at height zero qo decreases (more attraction), there 
is a threshold , which depends on p, such that the transition becomes first order for q^ < q^. In both 
cases the critical point q^ remains unaltered, see Fig. [T] 

In the bKPZ- class {O < p < 1) below gg, a new critical valué qc arises, this is shown in Fig.[7l For 
qo < 9o there is a phase coexistence región, in the sense that, dcpcnding on the initial conditions, the 
interface will be a moving or a bound one. For cxample, if a flat interface at a height far enough from 
the substrate is taken as initial condition, the interface will grow and not stay pinned to the substrate. 
Insidc the phase coexistence región, the bound phase is the stablc onc in the thcrmodynamic limit, 
i. e., the average time for the interface to detach from the wall grows exponentially with the system 

(2) 

size [26,31,32]. The critical behavior of this new transition, taking place at valué qc , will be addressed 
in Sec. 4. 

Considering the bKPZ equation ([6]) instcad of the RSOSW model, the situation is the same as the 
one depicted in Fig. [Ti with the parameter a playing the role of q, the parameter b the role of qo and 
c with a fixed positive valué. Note however that a and q ov b and q^ do not have precisely the same 
meaning and their exact relationship is not known. 

In equilibrium wetting, at a = Oc, the transition that occurs when approaching bw (analogous to 
(7o), for b < bw, is known as critical wetting. Whereas the transition taking place, for b > bw, by 
approaching Oc, with a < üc, is known as complete wetting [7]. Therefore, critical (complete) wetting 
corresponds to 9 = ge (90 > 9o) and q^ [q) approaching q^ [qc) with q^ < q^ (q < qc) in the RSOSW 
model. As pointed out in [32], the transition taking place for qo < q^ by varying q at the new critical 

point qc > qc for the bKPZ- class is not a wetting transition but rather a dcpinning transition 
because there is no phase coexistence at criticality (phase coexistence corresponds to q ~ qc, where 
the interface vclocity is zero). 

At the tricritical point, qo — qo - new critical behavior is observed. We define the exponents associ- 
ated to it using the supcrscript t. For example, at qo = qo, we have 



(29) 



p-^-iq-qc)" ■ 

We also define the foUowing exponents associated to critical wetting, 

Po - ÍQo - 10)'^ > 



(30) 



(31) 



w 



{qo - qo) 



(32) 



where the above relations are valid at q = qc. 
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bEW and bKPZ+ cases 



bKPZ- case 



le 




Fig. 7 Representation of the situations generated by the an additional attractive forcé between substrate and 
particles, for a fixed valué of s. For the bKPZ+ and the bEW cases the transition goes from second order 
(full hne) to first order (doted hne) depending on the valué of qo. The critical point qc is not changed. For 
the bKPZ- case there a is phase coexistence (PC) región. Considering the bKPZ equation, the situation is the 
same with a playing the role of q and h the role of go. 



II 



q 
III 



IV 



q ^ I ^ I I q 

Fig. 8 Four different interface configurations forming a closed cycle and the respective transition rates. 



3 Exact results 



The solution of the RSOSW model, when it is on the bEW universality class, is presentad below. 
For p = 1 detailed balance holds and the model can be solved exactly for {q < 1) [16]. First we 
solve the model for qq = q, then we consider the more general case < q [26]. Concerning exact 
Solutions obtained with the height probability distribution in the stationary state in the bound phase, 
our presentation foUows [33]. Also, by using a method introduced in [34], we calcúlate the velocity of 
a free interface for the RSOS model, Icading to an exact calculation of the exponent (3^^ which, to our 
knowledge, is not calculated elsewhere. 



3.1 Transfer matrix formalism 



In a general dynamical system, detailed balance is fulfilled if, for every pair of microscopic a and cr', 
the probability currents cancel each other, i. e. 

PaWa^a' = Pa'Wa'^a, (33) 

where is the probability of being in the state a in the stationary state and w^^^r is the transition 
rate from a to a'. From Fig.[8]we see that, for the RSOSW model, it is satisfied only if 

Pl = q-'PlI = q-^PlII = q-'p-'PlV = P^'ñ, (34) 

which implies in p = 1 . 

The detailed balance condition with p = 1, gives the foUowing equation, 

P{{K})^qP{{h,~l}), (35) 
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where thc configurations {hi} ~ hi, .., and — 1} = hi, /i¿ — 1, .., are such that the 
RSOS constraint is obcycd. A simple ansatz arising from this condition is 

L 

P{{h,}) = Zl'l[q'^^, (36) 

¿=1 

where is the partition function. It is given by 

L 

Zl=T.I[i''^ (37) 
{h} t=i 

where the sum is over all configurations satisfying thc RSOS constraint. Wc can writc thc probabihty 
distribution (|36p in the form 

L 

P{h,,...,hL) = Zl'l[n,,h^^„ (38) 

i=l 

where Th^^hi+i is the transfcr matrix and wc are using pcriodic boundary conditions (/il+i — hi). Thc 
elements of the transfer matrix are given by 

Tk,k' = q^^+^'^/^iSk^k' + Sk,k'+i + 4,fc'-i), (39) 

where A; > O and fc' > 0. Note that thc transfer matrix is tridiagonal and this is a dircct consequence 
of the RSOS constraint. 

Within thc transfer matrix formahsm. thc partition function rcads 

= E ■■■T.^inMTh,M-Th,^uhLTh,M, ^Y.^t,M = Tr(T^), (40) 

hi IiL hi 

where thc operator Tr gives thc trace of thc matrix. Wc are intercsted in calculating thc density of 
sites at height k, it can be written in thc foUowing form, 

L 

Pk = Zl'^'^Sh^^WTh^M+i = ZL^^--^Tkji,^...Th^,k = Zj^^T^',,, (41) 

{h} i=l ¡12 Hl 

where the tcrm 5hi,k imposes the constraint of summing only over configurations with hi = k. In 
bra-kct notation thc same quantity becomes 

pu^Zl\k\T^\k), (42) 

where the vcctors |fc) form thc canonical basis in height space. 

In most of the foUowing calculations the limit L ^ oo will be taken: on it 

Zl « A'^ (43) 

and 

with A being thc máximum eigenvalue of the transfer matrix and 10) the corresponding eigenvector, 
i.e. 

T|0)-^|0). (45) 

After setting up the transfer matrix formalism wc proceed to calcúlate the critical exponents of the 
bEW universality class. 
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3.2 Calculation of the critical exponents 



Based on dimensional analysis, (|2ip and (p4)) wc observe that the conditional corrclation function, 



is expected to follow the relation 



c(0 = 77 — (46) 

c{l) ^ r^^/''^, (47) 

vahd at the critical point. This quantity is equal to the sum of all possiblc paths, cach of which 
multiplied by it's respective weight, connecting two points at height zero with a distance l from each 
other. At the critical point all non-zero elements of the transfer matrix are equal to one, implying that 
all possible paths have the same weight. Therefore, for large enough l and at criticality, the number 
of all possible paths connecting two points at height zero at a distance l from each other is equivalent 
to the probability that a random walk starting at the origin will come back to the origin for the first 
time at time l. It is known that the probability distribution of the time that a random walker takes 
to return to the origin for the first time r is given by P{t) ~ T-3/2 [35]. From this foUows that, at 
criticality, 

c(0 - r3/2 (48) 

and, therefore, P/i'± = 3/2. We note that this kind of random-walk argument was successfuUy used in 
equilibrium wetting [36] . 



We now consider the problcm in the thcrmodynamic limit, L ^ oo. From cquations (|39p and (|45p 
foUows that 

g'=(g-i/20fc_i + 0fe + g+i/2^fe+i) = \(t>k, (49) 

whcre (j)k is a component of the vector ji^). In ordcr to solvc this cquation we take the continuum limit, 
which should be valid when we are cióse enough to criticality. On this limit: k h, (pk 4'{h) and 
equation ([49|) becomes 

-jr^(p{h) + (3 - A)(p{h) - 3eh(f>{h) = O, (50) 
where e = í — q. The solution of it is 

where Ai{x) is the Airy function. From the condition that (f){h) has to vanish for h < O foUows that 
yl = 3 and 

(j){h) ^ Ai{3^/^e^/^h). (52) 

With this explicit form of (^(/i) the calculation of the mean height and the interface width is straight- 
forward, they are given by 



{h)=A-^ 4>{hfh^^''^ w^Ja-^ 0(/^)2[/l-(/l)]2^el/^ (53) 



where A = (f){h)'^dh. Since 0(0) = O, in order to calcúlate the density of sites at zero height in the 
continuum limit we have to calcúlate p{h) at some small fixed height h = 5. Considering that 5 is much 
bigger than e and is small enough such that (j){5) « (5(/>'(0) is a good approximation, we have 

p(O)^A-i(0'(O)2), (54) 

which gives 

p(0) ~ e. (55) 

These exact rcsults are all for the stationary state, the present method does not allow us to calcúlate 
time-dependent quantities. In order to obtain them we resort to Monte Cario simulations [33] or 
numerical integration of a set of equations obtained with the supposition that the time-dependent 
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probability distribution is pair-factorized, which seems to be the case (see Sec.[Ü]). Both methods give 
thc foUowing rcsults, 

{h) ^ t^/\ w ~ ¿1/^ (56) 

and 

Po-t-'^'- (57) 

We havc obtaincd thc critical cxponents of thc bEW univcrsahty class for d = 1, they are: /? = 1, 
iy± — 2/3 and í^y = 4/3. The above calculations are not possible for the bKPZ case {p ^ 1), where we 
have to use numerical simulations to obtain the critical cxponents. Next wc consider thc ncw sccnario 
generated by an attractive substrate for the bEW class. 



3.3 Thc case qo ^ q 



With qo < q, below a certain valué of the deposition rate at height zero q^ thc phase transition 
becomes first-order. We now calcúlate the valué (jg, show that the transition is first-order for qo < qQ 
and calcúlate the cxponents associated to the tricritical point. 

The only elements of the transfer matrix that are changed, when qo ^ q, are Tq^q = q/qo, 7o,i = 
(g/qo)^^^ and Ti^o = [q/qoY^'^- The new transfer matrix can be written in the form 

n,k' = 9('+''^/'(9/9o)(''='°+'^''°^/'(4,fc' + + (58) 

where fc > O and k' > 0. Using (j38p . one can verify that the probability distribution obtaincd with this 
transfer matrix satisfies detailed balance. 

For q > 1, the probability distribution is clcarly not normalizable, independently of the valué of qo. 
This shows that the critical point is unchanged with the attractive forcé between the substrate and 
thc absorbed particles. For g = 1, wc make the assumption that 

= x'^ with fc > 1, (59) 

where a; < 1. Using the transfer matrix (|58p . thc set of equations obtaincd by applying (|59p to equation 
is 

9o ^^^<^o + X + x'^ ~ Ax, 



It has thc solution 



2(1 - qo) 2 qo 



x-^^ + l + x = A. (60) 



v/1 + 2^0 - Sq'é 1 



Which, with equation (¡^^ gives 



^ l + go-69^ + Vl + 2go-3go 
^° 2 + 49o-6qo ^ ' 

From this result we see that po = O at q^ = 2/3, po > O íor qo < qo and for qo > q*) thc assumption 
(|59p is not valid. Sincc po is finite at thc critical point q = 1 for qo < q^n thc transition is first-order. 



The cxponents related to critical wetting can be obtaincd cxactly as foUows. From equation (p^ . 
with qo below and near enough = 2/3, we have 

Po~(go-9o)\ (63) 
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giving Z?*^^' = 1 for thc bEW class. The dcnsity of particles at any hcight h can also be calculated, from 
(¡44¡) and jeT]), as foUows 

P'' = \ — TTTi 2\- (64) 

qo + x^/(l - x^) 

The above equation with (/i) = ^ kpi^ and vj^ ~ ^ {^)YPk, gives 

(/i) ~ ((Jo - go)"\ w-'(9o-9o)"\ (65) 

thcrcfore C^^' = 1- 

Thc critical cxponcnts defined at thc tricritical and it's vicinity, with q < Qc and qo — Qq, were 
obtained numerically in [33]. Thc resuhs are m agreemcnt with: 

(9c-g)'/', (66) 

{h)^{q,-q)-'/' w^{q,,-q)-''\ (67) 

Po ^ L-^'\ (68) 
whcrc thc first two relations are vahd ncar and bclow criticahty and thc third at thc critical point. The 
ofF-critical results can be confirmed by numcrical diagonahzation of the transfer matrix, whereas the 
finite-size result can be reproduced by evaluating of the product of L transfer matrices. Time-dependent 
resuhs, as in thc case of complete wctting, are obtained by MC simulations [33] or numcrical integration 
of thc set of equations, obtained with the supposition that the time-dependent probability distribution 
factorizes (sec Sec. 5), thcy are: 

{h) ^ (69) 

P^^t-^l^. (70) 
These results give the exponents /?* — 1/3, v^y^ ~ 2/3 and = 4/3, different from thc bEW critical 
exponents. Nevertheless, the exponents z — 2 and 7=1/4 are still the same. 



3.4 Exact calculation of the velocity of a free interface 



A method to obtain exact results, in the long time limit, for the RSOS model without the wall was 
developed by Neergaard and den Nijs [34]. Here we use this method to calcúlate the velocity of a free 
interface and, consequently, obtain the exponent . The calculations presented below are summarized 
in the appendix of [37] for a slightly different model. 

The variables ai = /li+i — hi can take only the valúes —1, 0, 1 becausc of the RSOS condition. With 
them the free interface problem can be mapped onto a problem of particles jumping in a lattice with 
the following rules: 



00 ^ 


> + - 


with rate q 


0+ 


+0 


with rate q 


-0- 


^ - 


with rate q 


-+ 


00 


with rate q 


00 ^ 


— h 


with rate -p 


0- 


-0 


with rate 1 


+0- 


^ + 


with rate 1 


+- - 


00 


with rate 1. 



Since the initial condition is a flat interface and all possible transitions conserve the total charge. the 
number of positive charges {a i = 1) in the system is equal to the number of negative charges (cr¿ = — 1). 
In [34] it was shown (in a more general context) that in the región of the phase diagram where 

2-9 
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the ansatz 

= , (73) 

where P{N) is the probability of having a configuration witli N positive charges in the hmit t —f oo 
and Zl is a normaUzation constant, is in agrccmcnt with the dynamical rules (j7ip . The normahzation 
condition gives 

^^=J^^NimiL-2N)\[2^) ■ ^^^^ 

In the thermodynamic limit, L — > oo, the sum in the partition ñmction is dominated by the term 
N = L/{2 + l\/q/ (2q — 1)), giving the foUowing density of positive (or negative) charges: 

-o(l/L). (75) 



2 + V(2g-l)/g 
The interface velocity is given by 

V^iq- p){00) + q{0+) - (+0) - (-0) + q{-0) + q{-+) - {+-). (76) 

Because of (00) = (1 - 2p)^ + o{l/L), (+0) = (0+) = (-0) = (0-) = (1 - 2p)p + o{l/L) and 
(-+) = (+-) =p^ + o(l/L). Henee, 

v = {q-l)2p + o{l/L), (77) 

where p is given by (j75|) . 



We obtained the asymptotic velocity for the free interface case, which is equal to the asymptotic 
velocity of the RSOSW model in the moving phase. Even though the exponent /3„ is defined in the 
horizontal direction in the qy~p plañe we expect it to be the same in other directions, because numerical 
calculations of the exponent (3y for different valúes oip are all compatible with (3^ — 1 [16]. Therefore, 
we have computed the exponent (3v exactly. 



4 Numerical simulations 



Here we discuss some technical aspects of Monte Cario simulations of the RSOSW model and summarize 
numerical results for the bKPZ- universality class obtained in [38] . We then establish a scaling relation 
between the exponents, based on heuristic arguments and in agreement with the numerical results, 
which shows that with the introduction of the wall just one new independent exponent arises. We 
explain how to intégrate the bKPZ cquation numerically in an eíRcient way and demónstrate numerical 
results bascd on this method. The attractive wall case and extensions of the problcm are also discussed. 



4.1 The exponents of bKPZ- universality class 



In order to calcúlate the critical exponents numerically by simulating a latticc model one can use 
ofF-critical, finite-size and time-dependent simulations. With off-critical simulations it is possible to 
calcúlate the exponents ¡3 and C, while the exponents 9 and can be obtained obtained with time- 
dependent and finite-size simulations, respectively. 

A technical problem with the off-critical simulations is the EW-KPZ crossovcr, which manifests 
itself as foUows. For a modérate simulation effort, it gives the impression that the critical exponents (3 
and C (measured for a fixed valué of p) depended continuously on p, varying from the EW exponent 
/3 = 1 and C = 1/3, when p is ncar 1, to larger valúes as p gets smaller, where we are considering 
the bKPZ- class (O < p < 1). But this is not the case, when making measures near the equilibrium 
point one has to access regions closer to criticality in order to get appropriates valué for (3 and C, they 
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Fig. 9 Off-critical simulations. Density of sites with height zero po (left) and the interface width w (right) 
as functions of the distance from the critical point qc — 0.4295(1), with p = 0.001 and L = 4096. Data taken 
from [38]. 

first look like bEW exponents and then they crossover slowly to the bKPZ- valúes. Because of this 
crossover, simulations with small valúes of p provide better results, that is why in [38] the simulations 
were done at p = 0.001. The off-critical simulations are presentcd in Fig. [9] and from it results 

/?= 1.67(5), C^í^-ia^ 0.41(5). (78) 

The EW-KPZ crossover also takes place when considering the bKPZ+ case {p > 1), with the difference 
that in this case, in order to get better numerical results, one has to use large valúes of p. 

As for finite-size simulations, what is appropriate for the RSOSW model is to measure how the 
critical point varies with the system size, because, as observed in [38], this variation is very pronounced. 
From relation (f24|) it is expected that 

gc(oo)-gc(i)~i"'/"", (79) 

where qdL) is the critical point for a system of size L and qdoo) is the cxtrapolated valué. With the 
use of this relation, Fig. [10] gives 

^ 1.00(3). (80) 

From the scaling relation ([29]) and the above numerical result we get the exponent ( = 0.5(1), in 
agreement with the valué coming from off-critical simulations 0.41(5). 

With time-dependent simulations the exponent 9 is obtained by plotting uq as a function of time 
at criticality. From Fig. [TOl 

e = pivw = 1.15(3). (81) 

From the scaling relations (P7¡) . (P5¡) . the finite-size result v^^ = 1.00(3) and z = 3/2 one can see that 
9 = 1.15(3) is in agreement with the off-critical result (3 = 1.67(5). 

Results coming form finite-size and time-dependent simulations are certainly more reliable than 
when extracted form off-critical simulations. Nevertheless, off-critical simulations are important to 
confirm scaling relations. Next we propose a scaling relation and show that it is in agreement with the 
numerical results presented here. 



4.2 Scaling picture and differences between the bKPZ universality classes 

We now show, with an heuristic argument, that with the introduction of the wall just one new inde- 
pendent exponent arises. We suppose that the velocity of the interface in the growing phase foUows 
V q — Qc, which is in agreement with numerical results and an exact result presented in the last sec- 
tion. Since it is the time derivative of the mean height, which has the same dimensión as the interface 
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L 



Fig. 10 Left: DifFerence A between the finite-size critical point qc{L) and the extrapolated critical point 
qdpo) — 0.4295(3) as a function of L for p — 0.001. Right: po as a function of í at the critical point Qc — 
0.4295(1), with p = 0.001, L = 8192 and 800 realizations. Data taken from [38]. 



case 


z 




H 









DP 


1.58 


1.10 


1.73 





0.159 


0.276 


bKPZ- 


3/2 


1 


3/2 


1/2 


1.184(10) 


1.776(15) 


bEW 


2 


2/3 


4/3 


1/3 


3/4 


1 


bKPZ+ 


3/2 


1 


3/2 


1/2 


0.228(5) 


0.342(8) 



Table 1 List of the critical exponents. Most of DP exponents in the first line come from [45], the exceptions are 
OL and C, that come from [40]. The bEW exponents come from the exact results presented in Sec. 3. The bKPZ 
exponents come the numerical results for the exponent Q obtained in [29] and the scaling relations presented 
here. Table taken from [38] 



width w, and has the dimensión of time, we expect that w ^ {q — qc)^\\ - With this, (¡22|) . p4|) and 
([27|) we have the scahng relation 

= (82) 

z ~ a 

In this argument we combincd relations valid above criticahty (w ~ g — Qc) with relations vahd below. 
However, one can see in table [TJ where the exponents of the bKPZ and bEW universality classes are 
displayed, that the seaUng relation ([82]) leads to exponents in agreement with the numerical (bKPZ 
classes) and exact (bEW class) results obtained here. A possible reason for that is: without a wall for 
q < Qc the interface would have a negative velocity, going linearly to zero as a function of the distance 
from the critical linc. It secms that, with the prcscnce of the wall, quantitics with the dimensión of 
the mean height divided by quantities with dimensión of time still go to zero linearly as criticality is 
approached. The scaling relation ([82]) . was obtained for the first time in [15], for all dimensions, using 
a diffcrcnt argument. 

With the scaling relation (j82|) and z = a/"f = v^/v^^ we have that the critical exponents i^y and i'± 
are determined by the scaling exponents. Thereforc, thcrc is only one independent (from the scaling 
exponents) critical exponent left, which is (3. 

All the above discussion should also be valid to the bKPZ+ universality class, with the diffcrence 
that the new critical exponent /? (or 6 = í3/i^\\) has a different valué. The best way to obtain the new 
critical exponent numerically is to pcrform time-dependent simulations with the SSW model, sincc in 
this model the critical point is known exactly. This was done in [29] , and the estimated exponents are 
summarized in table [TJ 

The exponent /3 is bigger (smaller) than one for the bKPZ- (bKPZ+) universality class, this tell us 
that a typical interface configuration of the bKPZ- universality class, at criticality, is characterized by 
a smaller number of contact points with the substrate, in comparison to a bKPZ+ typical interface. In 
Fig. [TT]we show typical interface configurations, obtained in [29] with the SSW model at p = O (A > 0) 
and p = 1 (A < 0). For the bKPZ- case, a smaller number of contact points and larger detached regions 
(distance between two contact points) are observed. 
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X>0 



Fig. 11 Upper row: typical spatio-temporal configuration of contact points with the substrate. Lower row: 
final interface configuration ai t — 5000. The figures were obtained with numerical simulations of the SSW 
model atp = 1 (A < 0) and p = O (A > 0). The number of contact points in the A < O case is clearly smaUer. 
Figure taken from [29]. 



Remarkably, the distance between the contact points difíers from the correlation length. More 
specifically, since the average distance of contact points l{t) is proportional to p^^it) and vj_ = 1, we 
have that l{t) > £,±{t) {l{t) < £,±{t)) for the bKPZ- (bKPZ+) class. The interplay of these two length 
scales leads to othcr difFcrcnccs between the bKPZ+ and bKPZ- universality classes, they are related 
to the distribution function of the distances between contact points and to the first depinning time 
probability distribution [29]. 



4.3 Wetting as a Contact Process with an external field 



Let US now turn to the special case p ~ 0. Herc the RSOSW model exhibits a very different dynamics, 
which was introduced in [39] and further analyzed in [37,40]. This is a very particular case of the 
model, because, with p = O, once a layer is completely fiUed an evaporation on it becomes impossible. 
Henee, since the initial condition is a flat interface at height zero, the presence of the wall makes no 
difference. 

It was shown in [39] that the phase transition at p = O pertains to the dircct percolation (DP) 
universality class [41.42], which is the most prominent universality class of nonequilibrium phase 
transitions into an absorbing state [43,44]. An explanation for this comes from the fact that the 
RSOSW model at p = O can be related to the contact process (CP) [45]. The CP is a well-known 
model in the DP universality class, it can viewed as a simple, and also simplistic, model for the 
propagation of a disease, where each site can be in two states: empty (healthy) or occupied by a 
particle (infected). Partióles can créate other partióles in empty first neighbors sites (propagation of 
the disease) or die spontaneously (cure) . If the system has no partióles (no sick individuáis) it is in the 
absorbing statc and the dynamics ceases. An non-exact map between the CP and the RSOS at p = O 
model was proposed in [39] where sites with zcro height are related to the infected individuáis in the 
contact process. Henee, once the first layer is filled it is analogous to enter the absorbing state in CP, 
since no partióles in this layer can be evaporated anymore and the height zero becomes inaccessiblc. 

In [38] it was shown that the región in the phase diagram with O < p < 1 can be interpreted as a DP 
process with an external field that dcstroys the transition. In the CP, an external ficld is introduced 
by allowing creation of partióles at a certain rate. With it, there is no absorbing state anymore and 
the transition is lost. In the same way the non-zero p destabilizes the absorbing state in the RSOSW 
model, since even after the first layer is completely filled an evaporation on it is still possible. With this 
interpretation the curvature of the phase transition line, Fig[3l is predicted and a crossover exponent 
from the DP to the bKPZ- universality class calculated (sce [38] for details). 
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Fig. 12 Numerical integration of equation (|86|l . Left: time-dependent simulations with the order parameter 
n = (m) as a function of time, at the critical point = 0.39025(25) for L = 2048, giving 9 = 0.235(10). Right: 
finite-size simulations with the saturation valué of n as a function of the system size L at the critical point, 
giving = 0.33(2). The valúes of the other parameters are b = 1, c = O, ct = 0.1, D = 1, s = A and dt = 0.1. 



4.4 Numerical integration of the bKPZ equation 



The MN equations can be integrated numerically very easily by using a method introduccd in [46] to 
intégrate Langevin equations with non-additive noise. The method consists of integrating the dcter- 
ministic and stochastic parts of a discrete versión of the Langevin equation separately in each time 
step. 

We now explain this method for the MNl equation in the case of a one-dimensional substrate. To 
this end we have to consider the spatially discrete versión of equation ([9|), which reads 

= -aui - - cnf"'^^ + cr[7ií+i + n¿„i - 27i¿] + UíQ, (83) 

where n{x,t) ni{t), x — !• i5x and we are using 5x = 1. The algorithm evolves as foUows. First a 
variable n* (í) is calculatcd from 

n*{t) = n,{t) + {-bn,{ty+^ - cn,{t)'^'+^ + a[n,+i{t) + n,^i{t) - 2n,{t)]}dt, (84) 

which corresponds to one step of the integration, using the Euler method, of the deterministic part 
(without the linear term) of equation ([55¡) . After that, n¿(í + dt) is obtained with 

ni{t + dt) = n*{t) exp{-adt + Voirj^), (85) 

which is the solution of the one- variable Langevin equation constituted only of the stochastic part and 
the linear part of equation (|83[) , where rji is a random number coming from a Gaussian distribution with 
zero mean and unitary variance. The reason to add the linear term in the sccond step of the update 
scheme is that the stochastic part can still be solved exactly with it. Note that the term — 2crn¿, 
included in the first step, could, instead, be included in the second step. After all the L variables are 
actualized according to the above scheme a step dt is completed. 

The numerical integration of the MN2 equation is more complicated because of the term (Vri)^/n 
[47]. In order to overeóme this problem, Al Hammal et al. [48] considered a non-order parameter 
Langevin equation. Such an equation is obtained with the variable m = 1/n that transforms the MN2 
equation (fTT|) into: 

d 

—m = am + brn^"'' + cm^~^* + ctV^to + mC,. (86) 
dt 

We integrated this equation numerically using the algorithm explained above. We performed time- 
dependent and finite-size simulations, obtaining 9 ~ 0.235(10) and /3/í^_l = 0.33(2) (see Fig. [T2|). The 
results are in agreement with the results obtained in [48] and the results presented in table [1] We note 
that changing the valué of s does not change the critical behavior [48] (in Fig.[T2]we used s = 4). 
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Fig. 13 Simulations of the RSOSW model with an attractive forcé between the substrate and the interface. 
Left: off-critical simulations giving (3 — 0.25(2). Right: time-dependent simulations giving 6 = 0.159(10). Both 
exponents are in agreement with the DP exponents. The simulations were done at p = 0.001 and qo = 0.39, 

where Qc = 0.4295(1) and g?' = 0.4377(1). 



4.5 Numerical results for the attractive substrate case 

As discussed in Sec. 2, in the bKPZ- case with an attractive forcé between the substrate and the 
interface there is a phase coexistence región if the attraction is strong cnough and a depinning transition 
takes place at a new critical point (see Fig. [71). 

This depinning transition was first obscrvcd to be first-ordcr [26,31]. Latcr Muñoz and Pastor- 
Satorras [24], with numerical intcgration of the MNl cquation, obtaincd exponents of the DP univer- 
sality class. In Fig. [13] wc prcsent ofF-critical and time-dependent simulations for the RSOSW that 
agree with the findings from [24]. Wc uscd the valúes p = 0.001, whcre qc = 0.4295(1), and qq = 0.39 

(2.) 

(which is smaller than í/q). Wc obtaincd the new critical point at q^' = 0.4377(1), /3 = 0.25(2) from 
off-critical simulations and 9 ~ 0.159(5) from time-dependent simulations, both are in agreement wit 
the DP exponents (see table[T]). We note that DP exponents were also obtaincd for other microscopio 
modcls [27,49]. Therefore, the obscrvation of a first-ordcr phasc transition sccms to be a transient 
effcct. 

Hinrichscn [50] argües about a possible connection between the pair contact process with diffusion 
(PCPD) and noncquilibrium wetting, pointing that the MNl equation is equivalent to the Langevin 
equation for the PCPD [51]. There has been a long discussion about the universality class of the 
PCPD [52], and it is still not clear if the PCPD model is in the DP universality class or if it is in 
a new universality class of modcls with an absorbing state. The result obtained in [24] would be in 
agreement with the PCPD bcing in the DP universality class, however, the continuum approach may 
be inadequate for PCPD [50]. 

Although the DP exponents were obtaincd in Fig. [T31 the transition. that takes place with the 
parameters used there, is different from the transition at p = 0. To see this one can consider the survival 
probability Ps (t) , which is the probability that the system will not enter the absorbing state until time 
t. At p = O the survival probability can be less than one because once po = O it will remain zero in the 
dynamics that foUows. It decays to zero, in a power-law (exponential) way at (above) criticality and it 
reaches a saturation valué after some transient below criticality. For p ^ O the survival probability is 
always onc (in the moving and in the bound phases) because even if the first-layer is completcly fiUed 
evaporation events on it are still possible. The Langevin equation related the RSOSW model at p = O 
is the DP Langevin cquation [26], where the noise tcrm is multiplied by the square root of the field, 
and it cannot be obtaincd from the MNl equation with a simple transformation of variables. 

The tricritical point of the bKPZ- universality class was analyzed, with numerical intcgration of 
the MNl equation and MC simulations of the SSW model, by Romera et al. [53]. They obtained 
7 = 0.35(2), z = 1.4(1), /3(2) = 1.50(9), C^^^ = 0.9(1), 6»* = 0.49(2) and = 2.0(2). As is the case of 

the bEW class, the critical exponents /3*, ly]^ and i'y diffcr from the critical exponents of the bKPZ- 
class while the scaling exponents are in agreement with the KPZ valúes 7 = 1/3 and z = 3/2. As far 
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as we know, the critical behavior of the tricritical point of the bKPZ+ universality class was not yet 
analyzcd. 

Wc point out that the wetting transition can be envisagcd as a contact process with long-range 
interactions, considering the contact points with the substrate as active sites and the detached part 
between two contact points as generating an efFective long range interaction between them [54-56] . In 
the so-callcd cr— process [54], dcpcnding on the paramcter controUing the long-range interactions the 
phase transition may be DP or first-order. 



4.6 Extensions of the problem 

Hitherto only short range interactions between the wall and the absorbed particles were considered. 
In order to study long range interactions one can take the potential 

^('^) = ¿ + ¿' 

where the parameters b and c have the same ñinction as they have in the short-range interaction 
potential (¡S]) and l > s. Hammal et al. [57] studied complete wetting for the long range interactions 
case with numerical integration of the associated multiplicative noise equation, coming form the Cole- 
Hopf transformation of the bKPZ equation ([6|) with the potential (|87p . and power counting arguments. 
The results are: for s > 1 the critical behavior is the same as the one obtained with short range 
interactions for the bKPZ- and bKPZ+ universality classes; for s < 1 mean field (in the sense that 
the noise is irrelevant) critical behavior holds. These results for long range interactions were also 
confirmed with a microscopio model in [57]. Another microscopio model with long range interactions 
was previously studied in [49], nevertheless, as pointed out in [19], it does not display a mechanism 
that produces surface tensión (like the RSOS constraint for example) and. therefore, it is not clear if 
it corresponds to the bKPZ equation with the potential ([87|) . 

One relevant point in considering long range interactions is that it may play an important role in 
possible experimental realizations of nonequilibrium wetting. Another generalization of the problem, 
that can be central in an experiment, is to consider diffusion of single atoms in the interface. This 
generalization was studied in [58] using the RSOSW model with diffusion. In this new versión of the 
model, another possible transition, taking place with rate D, is /i^ — > /i¿ — 1 and hj — > hj + 1, where j 
is one of the nearest neighbors of i, before the transition /i¿ > hi±i and the final configuration respects 
the RSOS condition. What was found in [58] is that, for p 7^ 1 diffusion can shift the critical line but 
the critical behavior is still the same and at p = 1 the critical behavior and the critical point do not 
chango. 

Another extensión, which is not yet studied, is nonequilibrium wetting with disorder. In equilibrium 
wetting, the random substrato and also random bulk cases were studied, and new physical propcrties 
are born from these situations [7]. An opon problem, in the equilibrium case, is what happens to critical 
wetting when a random substrate is considered [59-61]. 



5 Mean Field Approximations for microscopic models 

In this scction we present mean field approximations for the RSOSW and SSW models. We first 
write down the master equation for surface growth models and then apply simple and pair mean 
field approaches in order to solve it. The first mean field approximation, for microscopic models for 
nonequilibrium wetting, was introduced in [33] . The presentation here follows the works of Ginelli and 
Hinrichsen [28] for the SSW model and Barato and de Oliveira [63] for the RSOSW model. We also 
obtain some new results concerning pair mean field approach for the RSOSW model. 
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5.1 Master equation 



A stochastic procesa with continuous time evolves according to the master equation (scc [62]), such 
equation reads 



dt 

a' 



where Pa(t) is the probability of being in a configuration a at time í and w^'^a is the transition rate 
from the configuration cr' to the configuration a. In the case of surface growth models with deposition 
and evaporation rates depending only on height of the target site and it's nearest neighbors, the master 
equation becomes 

d ^ 

—P{hi,h2,...,t) = ^^{w„{hi^i,hi -n,hi + l)F(/ii, /i¿ -n,...,t) 

n i—1 

-«;„(/ií_i, /ij, hi+i)P{hi, hi, í)}, (89) 

where w„(/ii_i, hi—n, ft.¿+i) is the rate for a transition from (..., /i¿_i, hi—n, /i¿+i, ...) to (..., hi, /i¿+i, ...), 
L is the System size and we are considering a one-dimensional system with periodic boundary condi- 
tions. 

Denoting w„(/i¿_i ~ k,hi = l, hi^i = m) by Wn{k, l, m) we have that for the RSOSW model, where 
the height changes by only n = ±1, the rates are: 

w+{k, k, k) — w+{k + 1, k, k) — w+{k, k,k + 1) ~ w+{k + 1, k, fc + 1) = q, (90) 

w-{k,k + l,k) = w^{k,k + l,fc+ 1) = w-{k + 1, A: + 1, A:) = 1, (91) 

and 

w_(fc + l,fc+l,/s+l) =p. (92) 
On the other hand, for the SSW model they are non-zero only for n — ±2 and given by 

w+2{k,k^l,k) (93) 

and 

w-2{k,k + l,k)^l-p. (94) 

In the foUowing we denote P{k,l,m,t)wn{k,l,m) by Jn{k,l,m,t). The time evolution of the one- 
site probability distribution is obtained by summing over all heights but one in equation which 
gives 

^P(Z,Í) ^^{J„(fc,Z-n,m,í) - J„(fc,/,m,,í)}, (95) 

k,m n 

For the the pair mean field approach (see below) we also need the time evolution of the two-site 
probability distribution, it is given by 

^P(/í:, ^, í) = ^ '^{Jni'm, k- n,l,t) + J„(fc, l- n,m, t) 

m n 

-Jn{m,k,l,t) - Jn{k,l,m,t)}. (96) 



The scaling exponents of the KPZ and EW universality classes are related to the interface width w. 
One can make the stronger assumption that the one-site probability distribution of these universality 
classes, for an infinite system and in the long time limit, is given by 

p{Kt) = t-yO^\ (97) 
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whcrc f{x) is a scaling function. Sincc in thc calculations that foUows we always considcr an infinite 
system in the long time limit, when possible to solve the problem exactly, we use this ansatz for the 
one-site probabiHty distribution at the wetting transition and in the moving phase. 

The mean ficld approximations we apply to thc RSOSW and SSW models in the foUowing consist 
in approximating the probabiUty distribution P{hi, h2, ■■■,t) by a (simple or pair) factorized form so 
that the master equation bccomes tractablc. Thcy are apphcd to onc-dimensional models and expccted 
to capture some features of them, thcy are not expected to become valid above somc critical dimensión. 



5.2 The RSOSW model 



5.2.1 Simple mean field approximation 

In simple mean field, terms like P{k,l,m,t) are approximated by their factorized form, i.e., 

P{k, l, m, t) = P{k, t)P{l, t)P{m, t). (98) 

Note that this approximation does not take the RSOS constraint into account. For example, the 
probability P{k, fc + 2, fc, í) is zero in the original problem and in the simple mean field approach it is 
simple given by P{k, k + 2,k,t) = P{k, tfP{k + 2, t). 

With this approach the master equation for thc one site probability distribution (PS)) for the RSOSW 
model acquires thc following simplcr form, 

jPk = q{PLl - Pk) + (1 - 2(?)(F|Pfe+i - PLlPk) + 
(2 - q){PkP^+, - Pk-iPl) +p{P^+i - (1 - Skfl)PÍ), (99) 

where Pk denotes the one-site probability distribution, fc > O and P_i = O (in order to account for 
the presence of the wall at height zero). The factor (1 — Sk^), multiplying Pj^, comes from thc fact 
that evaporation at zero height is forbidden. Thc initial condition, corresponding to an initially flat 
interface, is P^ = S^ q. 

Below criticality we assume that Pk decays exponentially, 

Pk=Ax\ (100) 

whcrc A ~ í — X is a normalization constant and a; < 1. This ansatz is valid below thc critical linc 
where the interface is bounded, while at criticality x ^ l. By substituting (jlOOp in equation ([99]) we 
obtain 

- px^ + {q- 2)x^ + {2q - l)x + q = 0, (101) 



which, with X = í, gives the critical line 



9c = ^(p + 3). (102) 



With thc probability distribution (|100p wc have Pq ^ 1 — x, (h) = and w = From equations 
(|101|) and (|102|1 . one can verify that near thc critical linc 



2(ge - q) 
x^l-^^. 

Therefore, 



(103) 



Po^{qc-q)\ {h)^{q,-q)-\ w ^ {q, - q)-\ (104) 
giving P — 1 and C = 1- 
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In ordcr to solvc at and abovc thc critical line wc takc thc continuum limit, whcrc Pk{t) 
P{h,t) and h = kS. In this limit, to sccond ordcr in S, cquation ((99|) bcconies 



dtPih, t) = -5l2{q - qc)P'dhP + S'{2q+1 + 3p) [PidhPy + -P'diPj . (105) 

In the long time limit higher order terms are irrelevant even with á = 1. To see this one can carry 
out the calculation that foUows with general 5 and then verify it with the final result. Since we are 
interested in this limit wc procecd with the calculations setting 5 = 1. 

At the critical point the first term on thc right hand side of thc above equation vanishes, using the 
ansatz (|Ü7)) wc obtain consistcncy only if 7 = 1/4 and u = 0. The difFcrcntial equation for thc scaling 
function is 

f{x) + xf\x) + 2{2q + 1 + ?,p)[2¡{x)f{xf + f{xff"{x)] = O, (106) 



which has thc solution f{x) = \/l — 2(2g+i+3p) ' S'iving 



P{h,t) = Jt-y^ - ——^ — . (107) 

^ ^ Y 2í(2g+l + 3p) ^ ' 

Above the critical point only the first term in the right hand side of equation (jlOSp matters, the 
sccond is irrelevant for t —f 00. By foUowing the same proccdurc as in the prcvious case, we obtain 



7 = 1/3, w = O and 



f{x) + xf{x) ~ 36(g - q,)!{xfr{x) = O, (108) 



which has thc solution j{x) = ^x/12{q — qc), giving 



for q> qc 



Vcry surprising is the fact that with this very simple approximation we obtain the KPZ growth 
exponent 7 = 1/3 above the critical line. A problem with it is that we do not obtain the mean height 
growing linearly with time above criticality, i.e., v = O also for q > qc- Thc next step is to perform a 
improvcd approximation, that satisfies thc RSOS condition and givcs v > O in the moving phase. 

5.2.2 Pair mean field approximation 

In the pair mean field approximation P(k,l,m,t) = P{k,l,t)P{l,m,t)/ P{l,t). This mcans that thc 
probability distribution is approximatcd by a pair factorizcd form. Clearly, thc present approach sat- 
isfies the RSOS condition. 

For the pair mean field wc nced thc time cvolution of onc-site and two-sitc probability distributions. 
Because of the RSOS condition only Pk^k, Pk.k-i and Pk.k+i are non-zero. Thcrefore, wc havc four 
variables: the three two-sitc probability distributions and the one-site probability distribution Pk. They 
are not all indcpendcnt, onc obvious constraint is that Pk = Pk,k + Pk,k+i + Pk,k-i and the other is 
Pk,k+i = Pk+i,k - The second comes from the fact that the transition rates are symmetric and initially 
the number of up steps is equal to the numbcr o down stcps (ñat intcrface). Thcrefore, wc are left with 
two independent set of equations. 



Applying the pair mean field to equation (|96|1 we obtain 

'Vk-iiVk-i + Xk-i) Xk{xk + Vk) 



Pk-i Pk 



ykivk + Xk+i) XkVk-i 



Pk+i Pk 



x^ 



-2í7-f(l-4,o), (110) 
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Fig. 14 Comparison between the phase diagrams obtained with mean filed approximations, simple (SMF) and 
pair (PMF), and simulations (S). Figure taken from [63]. 



dt 



Vk 



vi 



vi 



P 



'-k+l 



k+1 



P, 



k+1 



(111) 



where Xk = Pk,k and yk = Pk+i,k = Pk,k+i- The term 1 — Sk,o multiplying 2p-^ in thc first equation 
comes from the fact that evaporation events are forbidden at height zero and the initial condition is 
Xk = 4,0 and í/fc = O for aU k. 

The above equations were integrated numerically in [63], the results are in agreement with [3=1, 
C = 1/3, the EW growth exponent 7 = 1/4 at the critical Une and the KPZ growth exponcnt 7 = 1/3 
abovc it. Also, with this improved approximation an interface growing hnearly with time for q > qc 
is observed. Thc phase diagrams, coming from mean ficld, are compared to the one obtained with 
simulations in Fig. (|15p . One can see that the agreement of pair mean field is much better when 
compared to simple mean ficld. 

In relation to the simple mean field the main improvements are: better agreement with simulations, 
when comparing the phase diagram and also observables (not showed here see [63]), and a interface 
growing with a non-zero velocity in the moving phase. The exponent C = 1/3 is different from the one 
obtained with simple mean ficld, which is ^ = 1, while the others are the same (including thc EW and 
KPZ growth cxponcnts at and abovc criticality respectively) . 

An intcrcsting fcaturc of thc pair mean ficld is that it bccomcs thc exact solution of the model 
when p = l. From equation (|35p . with thc pair mean ficld approach, foUows that 



q 



Xk Xk+l 



Pk Pk+i 



Vk = XkXk+i, 



which has thc solution 



Vk 



(112) 
(113) 

(114) 

This is equal to equation (j49p . with thc componcnts of thc eigenvector, associated to the máximum 
eigenvalue A, given by 0^ = \/Pk- 



^PkPk+i 

Equations (|112p and pi3p . with thc condition Pk = Xk + Vk + Vk-i, give 
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Fig. 15 Pair mean field results. Left: comparison of the one site probability distribution obtained from pair 
mean field (PMF) and from simulations (S) q = 0.99 and í = 500. Right: phase coexistence (PC) región at 
the (¡r X p plañe with go = 0.4 obtained fi-om numerical integration of the pair mean field equations. 



The pair mean field becomes the exact solution because, at p = 1 and q <\. the probability distri- 
bution is pair-factorized in the stationary sate. This comes from the fact that the stationary probability 
distribution can be written as a product of transfer matrices. A natural question that arises is whether 
the time-dependent probability distribution is also pair-factorized. If this is the case, the solution of the 
pair mean field equations (|110|) and would give the exact time-dependent probability distribu- 

tion. In order to check this, we compare the one-site probability distribution obtained from simulations 
and numerical integration of (jllOp and in Fig. [TS] Wc see that they are in agreement, suggesting 

that the time-dependent probability distribution is indeed pair-factorized. 

The case go 7^ ? can also be considered within the pair mean field approach. Equations (jllOp and 
(jllip . for fc = O and fc = 1, have to be modified to 

d „ a;o(a;o+yo) , ^- ^ 

— xq = — ^Qo evaporation part (115) 

d „ yaiva+xo) xi{xi+yi) . n^a\ 

—xi = ZQo 2q h evaporation part (Ho) 

dt Po Pi 

d '^n ^0 • /-1 1 ^\ 

-rVo ~ Qo-^^ — + evaporation part, (11 ') 

dt Pq 

where the evaporation part corresponds to the terms of equation (jllOp and (jllip that are not multiplied 
by q. By integrating these new equations we obtain: for q^ < the transition becomes first-order for 
p > 1 while for p < 1 a phase coexistence región (in the sense explained in Sec. [2]) is observed, and 
the transition at the additional line gi^^ > qc seems to be first-order. In Fig. \TE\ we show the phase 
coexistence región in the q y~ p plañe for go = 0.4 obtained with numerical integration of the above 
equations. The results presented in Fig. [Talare new. 



5.3 The SSW model 



For the SSW model, with the pair mean field approach, we are left just with one independent set 
of equations. This happens because, in comparison to the RSOSW model, there are the same two 
constraints and not four but three variables (P^, Pk,k+i and Pk^k-i)- We can define the height for 
an interface, satisfying the single-step constraint, in the following way With every pair of sites with 
heights given by /i¿, /i¿+i, we associate a variable /i¿ = min(/i¿, ft.¿+i), that can be viewed as defined in 
a point between the two sites. For example, the initial configuration, where /i¿ = O (/i,; = 1) if i is even 
(odd), corresponds to h'^ — O for all i. In this way the two-site probability distribution Pk^k+i can be 
viewed as a one-site probability distribution Pj, where j = k. In the calculations that foUows we use 
this new one-site probability distribution. 
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The mastcr cquation for the (ncw) one-sitc probability distribution within pair mean field for the 
SSW model reads 

^„ j pu pj 



Pj^p 



dt ' ^V^.-l+^.-2 Pj+Pj-1 

+ (1 - P) p '1" p ' ^ (118) 

whcrc the boundary conditions, to be specified bclow, depend on the velocity of the wall. 
The above equation in the continuum Umit reads 



dtP{K t) = -Sip - l/2)duP + ' 6^i-^ + 49,tP 

{dnPf ÁdnPfdlP [dlPf d^PdjP dÍP\ 
4P3 ^ 8P2 + 4P ^ 4P 8 y ^ ' 

where h = jS and we went until order á^. As we did in the simple mean ficld for the RSOSW model 
we set 5 = 1, becausc onc can show that higher order terms are irrelevant in the long time limit [28]. 

In order to solve equation (|119p we apply the ansatz ((97¡) to it. The first term produces a linear 
propagation of the intcrfacc, therefore 

v^p~l/2, (120) 

which is in agreement with equation (jl4p . At p = 1/2 the first two terms on the right hand side of 
(|119p vanish and a non-trivial equation is obtained only if 7 = 1/4, giving the EW growth exponent in 
the equilibrium case. For p ^ 1/2 the third term on the right hand side of (|119|) is irrelevant in the long 
time limit and a non-trivial equation is obtained only if 7 = 1/3, giving the KPZ growth exponent. 

We proceed, presenting the rcsulting diffcrcntial equations for the scaling function and thcir So- 
lutions for the cases p=l/2,p = and p = 1. In each case subtle boundary conditions have to be 
used, in order to account for the moving substrate. Our aim is to calcúlate the exponent 6, within this 
approximation. for the bEW and bKPZ universality classes. To avoid confusión we denote the scaling 
function / and the critical exponent 6, for each valué of p, by /p and dp. 

5.3.1 bEW case 

At p = 1/2 the differential equation for the scaling function fi/2{x) is 



(/1/2)'* - ^./l/2(./í/2)^/í/2 + (/l/2)^(^(/í/2)^ + fl/2fl/2^ 
+ (/l/2)'(^/í/2-^/í';2)+(/í/2)' 



O, (121) 



where x = ht . Integrating it we get 



^fl/2f ^ Jl¡2fl /2 
(/l/2)^ /1/2 



2xA/2 - + 2^^ - ri)^ = 0. (122) 



We have to solve the above equation with the appropriate boundary conditions. At p = 1/2 the 
substrate is fixed at height zero. therefore, the cvolution of the probability distribution Pj foUows 
equation (|118|) for j > O, while at zero height it foUows 



dt 2 2 \Pi +P: 

where the missing terms in the equation come from the facts that evaporation is forbidden at the 
substrate and Pj = O for j < 0. Now if we assume that the scaling function satisfies /'(O) < cxd we 
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Fig. 16 The scaling function f{z) for the bKPZ- (p — 1) and the bKPZ+ {p = 0) universality classes, 
obtained from numerical simulations (doted line) and mean field (MF) approach for the SSW model. Figure 
taken from [28]. 



have that Pj w t ^^^/(O) for j = 0, 1, 2. Substituting this in the last equation, we see that it is valid 
only if /(O) = O (the left side of the equation is proportional to and the right sido is proportional 

to 

The sohition of equation (|122p . satisfying the boundary condition /(O) = O, is [28] 



25/4 

/i/2(x) = ^a;2 exp(-xVV2) 



(124) 



which gives 0i/2 = 3/4. 



5.3.2 bKPZ 



For p ^ 1/2 the diffcrcntial equation we get is 



ifpY 



= 0, 



(125) 



where x ~ (h — vt)t . 

Obtaining the suitable boundary conditions, for the cases p = O and p = 1, that accounts for a 
moving wall, is more involved. With an argument similar to the one presented above it is possible to 
show that /i(0) = O and /o(0) ^ O [28]. With these boundary conditions the solutions of equation 
(|125p . that are physicaUy suitable, are [28]: 



P 



= 1 




if O < a; < So 
if xo < a; < oo 



and 



O 



/o(x) ex Az( 



(126) 



(127) 



where Ai{x) and Bi{x) are Airy functions and xq « 3.32426. The scaling function (|126p gives the 
exponent Oi = 4/3 for the bKPZ- case and (|127p 6*0 = 1/3 for the bKPZ+ case. They are different 
from of the exponents obtained from numerical simulations (see tablc [IJ . 

Although the mean field theory does not predict the correct valué of the critical exponent 9 in the 
nonequilibrium cases, it does predict differences between the bKPZ+ and bKPZ- universality classes. 
For the bKPZ- (bKPZ+) an exponent bigger (smaller) than the equilibrium exponent 9i/2 = 3/4 is 
obtained. Also, as is shown in Fig [161 the scaling functions, given in equations (|126p and (|127p . are 
qualitatively similar to the scaling functions obtained from numerical simulations. 
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With this mean ficld approach the one-site probability distribution for thc frec interface case can 
also be obtained, in this case the difFerential equations for the scahng functions have to be solved with 
diíFerent boundary conditions [28] . In [64-67] , several exact results were obtained for the free interface 
case with the polynuclcar growth model, which is in the KPZ universahty class. A very important 
qucstion is whcthcr, with thc mcthods used in [64-67], thc boundcd interface case can also be treated 
analytically, allowing onc to compute thc cxponcnt 9 for thc bKPZ univcrsality classes exactly. 



6 Nonequilibrium wetting in higher dimensions 



We now turn to nonequilibrium wetting in higher dimensions. First we present a mean field approach 
to the MN equations, which yields a qualitatively picture of what happens above the critical dimensión. 
Then we discuss the critical behavior that is expected in higher dimensions using power counting and 
rcnormalization group argumcnts. We note that the description wc prcscntcd for thc onc-dimcnsional 
case is much more complete and, therefore, if one wants a better undcrstanding on what foUows one 
should turn to the refercnces prcscntcd below. 



6.1 Mean ficld for the continuum model 



The mean field approximation to thc MN equations is done in the foUowing way [68] . The first step is 
to approach thc discrete laplacian by 

V^n, = ^(n^ - n,) « (n) - n„ (128) 

j 

where the sum runs ovcr all ncarcst neighbors and this approximation is expected to become valid 
above somc critical dimensión. The equation that results from it, in thc MNl case, is 

= ~an - bn'+^ - cn^'+^ + a{n - {n)) + nC,. (129) 

The associated Fokkcr-Planck equation [21,62] reads 

^P(n, (n),í) = -¿{[-(a - D/2)n - hn^+^ - cn"^^^ + a{n - (n))]P(n, (n), í)} 



1 



n''P{n,{n),t)]. (130) 



Solving the above equation in the stationary state we obtain 

P,(n, (n)) oc n2(a-.)/z.-i í ^b^s c,, _ 2a^\ 

\ US üs nü ) 

where Ps{ti, (n)) is the probability distribution in thc stationary satc. Finally the critical behavior of 
the order parametcr (n) is obtained with thc sclf-consistcncy condition 

, , f!^ nPs(n, {n))dn 

in) = ' ■ (132) 

Jo Ps(n,{n})dn 



For the MN2 case we consider the non-order parameter Langevin equation (|86p . with the approximation 
(|128p . The solution of the corresponding Fokkcr-Planck equation is 

P,(m, (m)) oc TO2(a-<x)/D-i í _ ^^^s _ ^^~2s _ Mm^A (^33) 
^ ^ \ Ds Ds mD ^ ' 
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With the mean field approach introduced we procced considcring complete wctting, therefore, we 
set c = O and b constant and positive. In the calculations bclow, for the MN2 case we foUow [48] and 
for the MNl case we foUow [69,70]. 

For the MN2 case, the self-consistency condition (I132|) with the probabüity distributfon (|133p leads 
to the foUowing equation, 

^ MR) ' ^'''^ 

where 

hiy) = J^°° dxx^+^^-^y^-^ exp - -^a;- - ^) . (135) 

We want to calcúlate the exponent /?, defined by {m) ^ («c — a)~^ {m is a non-order parameter). We 
can write Ik {y) as 
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(136) 



where z = At criticality (m) —s- cx), we are interested in the asymptotic form of Ikiv), with ?/ — > oo. 
In this limit; 



h{y) « (^-^ j r(-2a/i5 + 2a/ D - k) 



^^r{-2a/D + 2a/D-k + s), (137) 



where the term exp — -^[^■^¡j was expanded until first order in (z/y)^. The above asymptotic 
form for /fe (y) and equation (|134p lead to [48] 

(to) - (-i:>/2-a)-l/^ (138) 

which gives /3 = 1/s and «c = —D/2. 

For the MNl case the resulting cquation, coming from (|13ip and (|132p . is [69,70] 

Jo((?^)) 

where 

2b ^ 2ay\ 



^ . (140) 

Ds xD ^ ' 



roo 

Jk{y) = / exp 
It can be shown that, for y ^ O, [69, 70] 

leading to 

(n) ^ {D/2 - a)niax[i/s,D/2.] ^^42) 

where max indicates the máximum. This shows that for the MNl case we have a more complex 
critical bchavior: a weak-noisc rcgime, where D < 2a /s and /3 = 1/s, and a strong-noisc rcgime, 
where D > 2a /s and P = D/2a. With the asymptotic form (jl4ip one can also calcúlate the critical 
exponcnts relatcd to higher order moments {n^) = Jk{{n)) / Jk-i{{n)) . The strong-noise regime can be 
further divided into two regimes that are difFerent with respect to the critical behavior of higher order 
moments [70]. 

Critical wetting can be studied with the same kind of procedure, this was done by de los Santos et 
al. [71]. In critical wetting, for the MN2 case just one regime is found, with the critical exponent being 
independent of the noise strength, while for the MNl case weak and strong noise regimes are found. 
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Moreover. with thc prcscnt mean ficld approximation, phasc cocxistence for thc bKPZ- univcrsality 
class, when thc attraction betwcen thc substrato and thc mtcrfacc is strong cnough, is observed [32.72]. 

Wc point out that in contrast to the mcthod presented above a simpler approach would be consider 
thc one-variable case by taking out the Laplacian. The one-variable MNl equation was solved exactly 
in [73], it has a rich scahng behavior, nevertheless it does not display the strong noise regime. One 
interesting featurc of thc onc- variable approximation is that it allows onc to clarify cssential diffcrcnces 
betwcen thc MNl and thc DP Langcvin equations [74]. 



6.2 Scaling analysis 



Given a Langcvin equation wc can define thc partition function Z by summing over all configurations 
and realizations of noise that satisfy it. In the case of thc MNl equation ([Ü]) with the potential (|10p 
and c = O it reads 

Z ex / DnDC,P{C\8[^n + an + hjf - V^n - nc\ , (143) 



where 

P[C] oc exp(CV2L') (144) 

and J DnDC^ denotes a functional integration. With the introduction of a response field ñ one can 
intégrate out the noise [75,76], rcsulting in thc foUowing equation, 



Zcx I DnDñcxpl -S[n,ñ\], (145) 

with thc action S[n, ñ] given by 



S[n, ñ]= J éxdt 



— ñ^n^ — ñ( —n + an + biiF — V'^ít- 
2 V oí 



(146) 



where, for simplicity, wc set <t = 1. Preforming naive power counting in this action wc obtain 

[n\ + [ñ] = -d and [D\ = d ~ 2, (147) 

where [x\ represents the dimensión of the quantity x in units of Icngth. Therefore, the critical dimensión, 
above which thc noise becomes irrelevant, is de ~ 2. It is known that thc response field ñ scales as 
the survival probability [77]. In the case of thc bKPZ- univcrsality class the survival probability is 
always one, therefore the dimensión of the response field ñ is zero, giving [n\ — d. Power counting at 
thc critical dimensión de ^ 2 gives 

/3 = 1, iyj_ = 1/2, z/|| = 1. (148) 



Using standard methods [78] one can perform a perturbative expansión with the action (|146p and 
then obtain the renormalization group flow diagram for the MNl equation. The calculations can be 
found in [18], in thc foUowing wc discuss some important points of it. Since the term proportional to b 
goes to zero in thc moving phasc, thc bKPZ equation is cquivalcnt to thc KPZ equation in the moving 
phase, the fiow diagram is similar to the well-known KPZ one [79]: above the critical dimensión de ^ 2 
there is a wcak- noise attractive fixed point. where D = O, and a strong-noise rcpulsive fixed point. If 
the noise strength is smaller then a certain threshold the flow runs to the weak noise fixed point and 
the critical exponents are given by (|148|) if it is larger the flow runs to infinity and the critical behavior 
is not accessible through perturbation theory. Therefore the situation is similar to the one obtained 
with the mean field approximation. This was verified with numerical integration of the MNl equation 
for d = 3 in [68]. 

Another important point about the perturbative expansión is that changing s, the hardness of 
the wall, does not introduce new divergences and therefore it is not expected to affect the critical 



33 



bchavior [18]. This is in agrccmcnt with numcrical rcsults and diíferent from the mean field result that 
prcdicts an cxponcnt dcpcnding on s in the wcak noisc rcgime. 

The critical behavior at the critical dimensión d ~ 2 cannot be determined by perturbation the- 
ory. This is a very important case because in experimental situations the substrate is usuaUy two- 
dimensional. The critical exponents 9 (or /5) of the bKPZ universality classes in two dimensions are 
not known. 

The renormalization group flow equations can be derived in an alternative way without writing 
down the action, but dircct from the Langevin cquation. For the MNl equation this can be found 
in [20] and for the KPZ equation in [10]. 



7 Final remarks 

Just as the KPZ cquation represents a robust universality class of nonequilibrium growing frce inter- 
faccs, the bKPZ cquation is expcctcd to rcprcscnt a robust universality class of nonequilibrium growing 
interfaces in the presence of wall. While equilibrium wetting transitions can be studied with the bEW 
equation. nonequilibrium wetting transitions are dcscribed by the bKPZ cquation. Bclow wc point out 
what wc considcr the main opcn problems in nonequilibrium wetting. 

With the introduction of the wall critical exponents arise, just one of thcm is indepcndent whilc the 
othcrs can be determined from scaling relations and the KPZ scaling exponents. An important open 
problem in nonequilibrium wetting is the exact calculation of the exponent 9 for the bKPZ universality 
classes. While an exact solution for the free interface case, for a specific microscopic model in the KPZ 
universality class, is known [64-67] the boundcd case still remains without an exact solution. As we 
showed here within an mean field approximation for the SSW model the exponent 9 and the one-site 
probability distributions for the bKPZ universality classes can be determined analytically but thcy 
differ from the numerical rcsults obtained for the fuU model. 

As pointcd out, extensions of the problem that were already considered are long-range interactions 
between the substrate and the absorbed particlcs and the study of the RSOSW model with diffusion 
of particles. What was not yet studied is nonequilibrium wetting with disorder, which can be very 
relevant in a experimental situation. Also relevant in possible experimental realizations, is to consider 
nonequilibrium wetting in a two-dimensional substrate. The critical exponents of the bKPZ universality 
class in d = 2 were not yet determined. Performing Monte Cario simulations with microscopic models 
or numerical integration of the Langevin equation at d = 2 is a trivial task, the main problem is to 
find some approximativc method that can support the rcsults obtained with simulations. 

The main challenge, in what we defined here as nonequilibrium wetting, is to observe experimentally 
the critical behavior obtained theoretically. Any growing interface in the presence of a wall and under 
nonequilibrium conditions is, in principie, a candidate of an experimental realization of the bKPZ 
universality classes. FoUowing the discussion in [19] good candidates may come from crystal growth 
and synchronization transitions in extended onc-dimensional systems. 
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